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Volkov wave function: its orthonormality and completeness 


Enderalp Yakaboylu 

The present note aims to provide a clear and explicit derivation of the orthonormality 
condition, and the completeness property of the Volkov wave function. Thus, none of the 
results are new. 


I, VOLKOV WAVE FUNCTION 


The Dirac equation for a charged particle interacting with an arbitrary electromagnetic field is 
written as 

d\4f{t)) ^ ^ 

i =[a-{p- eA) + ecp + m/3] |i/f(0), (1) 

where we used the natural units such that fi = c = 1, e is the charge of the particle, a and /3 are 
the Dirac matrices chosen in the standard representation [HI, and the gauge potential is given by 
= (0, A) with the metric convention (+, -). For the position representation of the state 

vector the Dirac equation Q reads 

[iy + ieA^) - m] = 0, (2) 


with = (J3,J3a). 

If we multiple Eq. Q with the operator iy''Dy + m from the left side, we provide 


0 = [yOy + m] [y y - mj t) = 


2 , 


+ 


[Dy, DA + m 


(A(x, t), 


(3) 


with D^ = d^ + ieAfj being the covariant derivative, t) = and we used the commutation 

relation {y^, y} = Using the fact that [Dy, y ] = ieFy^, where y = dyA^ - d^Ay is the field 
strength tensor, Eq. Q can be further written as 


le 


d + lieA^d^ - e A + —yy y + m 


ifr(x, 0 = 0, 


(4) 


where we have used the Lorenz gauge condition dy^ = 0. We observe from Eq. Q that as the 
gamma matrices couple to the field strength tensor, the spinor part of the wave function is gauge 
invariant for an arbitrary electromagnetic field. 

Our point of interest is the solution of Eq. Q for a plane wave. A plane wave can be defined in 
terms of the following gauge potential A^ = A^{ri) with the phase rj = k^x^, where is the wave 
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vector fulfilling the condition = 0. For a plane wave the spin resolved term of Eq. (|^ can be 
given by 

—y^y^Fy^ = iey^y^dyA^ = iey^y^A^ky , (5) 

where the dot denotes derivative with respect to the gauge field’s argument rj. Now let us seek the 
solution of Eq. Q in the form of 

t) = exp , ( 6 ) 


where (+) and (-) label the particle and anti-particle solutions, respectively. The free particle/anti¬ 
particle spinor with the normalization = 6sy can be given by 


= 


jp^ + m 

2 p° 


p-o-<f>s 
\ + m J 


^p,l = 


jp^ + m 

2 p 0 


( P-o'Xs \ 
p^ + m 

Xs 


(7) 


where p^ = (p°,p) is the free particle/anti-particle’s four-momentum with the energy = 
> 0, (f>i = X 2 = (1 0)^, 4>2 = Xi = (0 1)^ ^re two component spinors, and cr 
are the Pauli matrices [[D. In the absence of the plane wave, reduces to the identity, and 

t) becomes the field-free solution of the Dirac equation. If we plug the ansatz (|^ into Eq. Q, 
the unknown function satisfies the following equation 


+ 2ieA^ [+iPpf~\r]) + 

+ iey^kpy'"AyfF{ri) + m^fF{r|)^^ uf) = 0 . ( 8 ) 


Before going further, we first observe that the Lorenz gauge condition for a plane wave can also be 
written as = d{Mkp)ldp, which implies A^kp equals to a constant, and in fact, without 

loss of generality, this constant can be set to zero. Using further the on-shell condition p^ = m^, 
Eq.[^ simplifies to 

\+2iifkpf-\p) + (±2eA^^pp - e^A^ -t tey^k^y^Ay) mJ,;] = 0, (9) 


and whose solution just yields 


= exp 



exp 



( 10 ) 


where ^ = y^kp, and 770 is the initial phase at which the gauge potential and the field vanish. We 
further notice that 


exp 



= \ +e 


2pk ’ 


( 11 ) 
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due to the fact that = 0. As a result, the wave function with the quantum numbers p and s reads 


t) = exp 


-i ±ffXa ± 


— J" drj' {±Ap - eA^/2j 


^4 \ (+) 

1 + e -' 


2pk 


p,S ’ 


( 12 ) 


which is called the Volkov wave function [l2l. 


II. ORTHONORMALITY OF VOLKOV WAVE FUNCTION 


In this section, following Ritus [j3l|, we explicitly show that the Volkov wave function obeys the 
usual orthonormality condition 

J (Px t) = 6{p' - p)6,'^s ■ (13) 

For a shorthand notation, we first rewrite the Volkov wave function in the following form 


ipflix, t) = exp [-/ (l + e 


Ipk 


,(±) 

^p,S ’ 


(14) 


with the phase 


= ffXf, + r dT]' (±Ap - eA^/2\ . (15) 

Jrjo 

Then, we set the scalar part of the gauge potential zero within the Lorenz gauge, and choose the 
propagation direction along the z-direction and the vector potential on the x - y plane, i.e.. 


A^iT]) = (0,AAt-z),Ay(t-z),0) . 
Thereby, the phase of the Volkov wave function is 


=P^t-PzZ 


-p,x+ f\ 

Jo 


, +2ep_y ■ A + 


di] 


(16) 


(17) 


/o 2p_ 

where for the sake of convenience we set rjo = 0, Pj^ = {px, Py), and p- = p^k^ = p^ - p^ > 0. The 
phase can also be rewritten as 

^t-z /- ^ „ A \2 


= A 


-PzZ-p^-X+ f 

Jo 


(p^ + eA) -pI t + z 


-p^x+ r 

Jo 


drj 


/ {Pi. + ^^) + 


2p. 


(18) 


As the integrand is independent from x and y, one can immediately evaluate the corresponding 
integrals in Eq.[T^ which are in fact Dirac delta functions. Therefore, the orthonormality reads 

(19) 


J d^x 4,yJ(x. t)py,(x. o = J d^x exp [±;(5y - sw)] V!3 


S(P‘ 


Z-Px)f 


dz exp 


+i 


i, y r '^ 
(p--p-)^ — (p_-p-) 


f 


2 \ 


drj' 


2p-p'- 


/j 


f;(±) 

^ p',s' ^ p,S ’ 
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Now, we define 


whieh implies 




, Ap± + eAf + m2 


+ 

and henee simplifies the integral to 


P-P- 

{pj^ + eAf + m2 
P-P- 


J (fx t)ij/f^^(x, t) = 6(p\ - pj 

X ^ ^ , - 

J-oo 1 + [p^ + eAj + niA Ip-p'_ 


^-oo 1 + + eA)^ + m^j / p-p'_ 

The spin resolved part can be written as 

and the intermediate square matrix can be decomposed as 


±l{p_-p-)—^ . 


1 f+) 

‘*v 4"“’ 


, ^ /«/W, ^ M 


1 , 1 / 2^2/(/-r^) , lipL -p-)y^ - {pL 

1 + ^— \e A ---+ 2eA,y --- 

P-P- \ 2 4 


+ p_)/ 


{p'_ - p-)y^ - (p'_ + p-)y 


+2eAyy‘ 


If we evaluate each corresponding term, we provide the following expressions 


u , , 

p',^y 

.PuP _ ,,3', 


(±) t i{p'--p-)y - ip'- + p-)y 

u , , y - 

P .V ' A 


4 

,3 _ ( w 


r-^ 

where the common function 



(27) 

i*] = gip',p)s s',s. 

(28) 

’f} = Pxgip',P)S s',s. 

(29) 

'2 = Pygip'^ P)Ss',s. 

(30) 


gip',P) 


ip'_ + m) ip_ + m) + p]_ 

4 y(m + p^)im + p'°)pV'° 
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As a result, we obtain 


= (l + [ip± + /P-P-)8(p\p)3s\s , 


which exactly cancels the corresponding denominator in Eq. (231, and therefore we find 

J = 6{p\ - pJ5 2 ) sip^ p')5s\s • 

The function g{p, p') reduces to 




for the same momentum. Finally, using the relation 

P- lp'_-p- 


2p^ 


S{p'y - Pz) . 


we show that 


f 


t) = 6{p' - . 


(32) 

(33) 

(34) 

(35) 

(36) 


III. COMPLETENESS OE VOLKOV WAVE EUNCTION 

In this section, we present a clear derivation of the completeness property of the Volkov wave 
function, explicitly we will show that 

C(x,x') = ^ r d^p{il/^pl{x,t)ifj^p^{x\t) + if/^~l{x,t)if/^~f{x',t)) = 6{x-x'). (37) 

.v=l 

Using the form of the Volkov wave function given in Eq. ( [T4] ), the completeness can be given by 

C(jc, jc') = J £p exp [-iS ^^\x) + iS ^+)(jc')] U^^XxW^pfix') 

+ J J^/7exp[/S(“^(x)-/S^“^(x')] C/J;j(x)f/J,7f(x'), (38) 


Ul,%\x)Uj,f(x') = 


1 ± e 




H(x) ] f ± 


m 


^ 2p_ ) 2p0 ^ 


y%{x'W 

2p'- 

yU{x')ly^ 
2p'- 


(39) 

(40) 


with 
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where we used the completeness of the free particle/anti-particle spinor in the second line. With the 
help of the phase ( [T8] ), Eq. ( [38] ) is rewritten as 

- eAf + 


C(jc 




(fpe 


—{z-z)-p^-{x-x 


-I 


2p. 


u^;}{x)u^;f{x') 


+ 




I 

e L 


■{z - z) -p^-{x-x')- f 

Jt-Z 


{Pj^ + eA)^ + 

2p. 


U^;}{x)U^-f{x')- (41) 


After the following transformation p ^ -p in the second integral, we provide 

{p^ - eAf + irP- 


.x') = fd^pe 

+ J d^pe 


Y(z - z') - P^ ■ (x - x') - J' 




~(z-z!)-P^-{x 


{p^ - eAf + rrd 

2p, 


Ufl{x)Uflf{x'), 


with /?+ = /)° + p^. Now, for the first integral we apply the coordinate transformation {p) 
and for the second integral {p) {p^, p+). In other words, for the first integral 


(42) 




Pz 


111 

m + P j^- p 


2p. 


and for the second integral 


Pz 


111 
m + p^- pX 


-2p, 


X OO 

dpz = I 

OO kJQ 

X OO /^O 

dPz = I 

OO ^0 


dp- 


+ p\+ pi 


2 pi 


(43) 


dp^ 


m^+pl+ pI 


2pi 


(44) 


Then, the completeness property becomes 


C(x,x') 


u^:Kx)u':f(x-) 


p , „2 , „2 \ -I 


2pi 


'-{z - z) -p^{x-x 


'>-r 


(Pj^ - eAf + rrd 

2p. 


X 


+ 


,2 , „2 , „2 \ -I 




2p 


'-{z-z')- p^-{x 


x')+ f 

Jt-^ 


(Pj^ - eAf + 

2p^ 


xf/i-;»t/«,V). 


(45) 


























If we further apply the transformation p+ -p_ in the seeond integral, we provide 
C(x,x') = 


a , \ —I 




2p£ 


-Yiz-z')-p^-{x-x 


nt-z' 

Jt-z 


(Pj^ - eA)^ + ird 

2p. 


X u^;j{x)u^;f{x') 


,2 , „2 , „2 \ -I 


f d^p^ f dp_ ( I e 
%J ^—00 




y (Z - Z) -P^ {X 


-x')- f ' 

Jt-Z 


(Pj^ - eAf' + trd 

2p- 


X u^:lxx)u^:iJ{x') 


p+=-p- 

At this level we observe that 


u^;](x)u^;]\x') = u[-lxx)u^sls\x') 


p,S V*'*'/ p^s 

therefore, the completeness just yields 
C(x, x') = J" d^p_^ exp {ip^ ■ (x - x' 

where we used = (m^ + p]_+ P-) /{2.p-), and 


p+=-p- 


(46) 


(47) 


p^oo 

1 dp- exp 

-l\Kip. - 

J —00 

\ p- 


^U^^}{x)U^;f{x'), (48) 


P- 


K\ 


z- z 


K2 


f 


{p^ - eAf + nd- 
2 ’ d,_, 2 

such that kiK 2 > 0. 

The spin resolved part ([47]), on the other hand, can be written as 


(49) 




p,S V'*'/ p^s 

if e 

+ m H- 


(50) 


^{A{z)t - A(z)l/t)(^ + m) - 2p-A{z) + 2pA{z)t - eW(z)4(z')] [ 7° • 


2p- 2p 

We notice that the above equation is in the form of 


^u^;}{x)u^;f{x') = iCo + — + ^ 


K, K. 


P- P- 


where 

7°+ 7^ 0 

Ao - - - -7 , 


1 


Ki = ^\m-p^-y + ^ [Aiz') - Aiz)] 7 ^( 7 ° + 7 ^) - eAiz') 
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(51) 

(52a) 

(52b) 


K 2 = ^\trA + p\ + e [Aiz) - Aiz)] p^-y + em {A{z) - A{z)\ - e2p^ ■ A{z) - e^A{z)A{z)\ ^ 7 **. 

(52c) 
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Finally, using the expression (see Appendix [TV]) 


X o 

Ot 


dp- 



./ K2\ 

exp 





6 {ki) n = 0 


0 


n = 1 5 


with 


6iKi) = 26iz-z'), 6iK2) 


6 {k 2 ) n = 2 

26(z - z') 

{p^ - eA) + nd 


we prove that the completeness property can be written as 




(53) 


(54) 


C{x, x') = J exp (ip^ ■ (x - x') ((/7^ - eAf + m^) - y^) + ^^^(/ + y^) 


d(jc - jc') 


r 


(55) 
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IV, APPENDIX 


We will evaluate the following integral 


X oo 

dxfn{x) = I 
00 kJ—& 


dx 



r ./ ^^2)1 

exp 

1 


(56) 


with Ki > 0, and K 2 > 0. Accordingly, we consider the corresponding complex integral along the 
contour shown in Fig{^ 

/C2' 


dz fniz) =(pdz 
c Jc 


exp 


I U<lZ 


0 . 


(57) 


where the second equality is due to that there is no pole inside the contour. The complex integral 
can be decomposed as 


I dzf„(z)+lini I , 

Jc„ Jcr 


dzfniz) = I„+ lim I dzfniz) + lim | dzfniz) = 0, 

p—*0'^ 


( 58 ) 
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then the integral /„ can be given by 


r idO I iK2 K2 

4 = - hm „ , . exp iKipcos(e) - Kipsm(9) - cos(e) -sin(0) 

\ r ide 


+ lim 

R—^oo 


r u 

Jo 


Idu / IR? Ki \ 

cos(0) - KiR sin(0) + cos(0) - ^ sm(e)j . (59) 


We first observe that when n = 1, the integral vanishes, Ii = 0. In the case of n = 0, the integral 
along the contour Cp vanishes, and the integral Iq becomes 


/o = lim I id 6 Re exp (-iKiR cos( 6 ) - kiR sm( 6 )) , 

S->oo _ L 


which vanishes as long as /Ci 4^ 0. In the case of K\ = 0, the integral yields 4 = limR_,oo 27?, which 
goes to infinity, as a result 4 = In a similar way, when n = 2, the integral along the contour 
Cr goes to zero, and we have 


4 = lim 


r ide 
Jo P 


e exp — cos(0) - — sin(0) 
\ P P 


which is zero for K 2 4 0, otherwise 4 = limp^o+ 2/p, i.e., 4 = 5 {k 2 ). To sum up, we conclude 


exp -i \Kix 


d(A:i) n = 0 


- 0 n = 1 • 


6 {k 2 ) n = 2 


Imiz) 



FIG. 1. The contour for ki > 0, and /C 2 > 0 is shown. In the opposite case, the contour can be chosen in the 
upper-half plane. 
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We note that in the ease of ki < 0, and K 2 < 0, whieh also satisfies kiK 2 > 0, one ean ehose 
the integration eontour in the upper-half plane of the eomplex plane, see for a detailed diseussion 
Ref. a. 
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